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Abstract 

We perform the Hamiltonian analysis for a nonprojectable Hofava model whose 
potential is composed of R and R 2 terms. We show that Dirac's algorithm for 
the preservation of the constraints can be done in a closed way, hence the algebra 
of constraints for this model is consistent. The model has an extra, odd, scalar 
mode whose decoupling limit can be seen in a linear-order perturbative analysis 
on weakly varying backgrounds. 

Although our results for this model point in favor of the consistency of the 
Hofava theory, the validity of the full nonprojectable theory still remains unan- 
swered. 



1 Introduction 



Since Hofava formulated his theory [TJ of a renormalizable quantum gravity based on 
anisotropic relativity, there has been a debate about its theoretical consistency. Several 
analyses have been performed on models having the anisotropic relativity symmetry 
based on Hof ava's proposal, some of them can be found in Refs. [21 (3J HI |5l El [7] . The 
main discussion has been focused in the number of physical degrees of freedom, the 
physical consistency of their propagation, the coupling/decoupling of extra modes at 
the IR and the structure of the constraints together with the closure of their algebra, 
which is the point that mainly concerns us in this papei0. 

One of the issues of the general discussion raises at the very definition of the action, 
since the gauge symmetry group of Hofava models allows to formulate them in two 
different ways: the projectable and nonprojectable versions. In the projectable case the 
Hamiltonian constraint is automatically a global constraint. As a consequence, such 
models have an extra scalar mode at each order in energy scale. Hence, there is not 
a clear way in which general relativity can be recovered in the projectable case, unless 
other ingredients are added to the theory [9]. Because of this we concentrate our study 
in the nonprojectable scheme. 

In a previous paper [7J we completed successfully the Hamiltonian analysis for the 
lowest-order effective action of the Hofava theory. This action, which we call the XR- 
model, is given by 



and it is the action one obtains if the potential of the original Hofava theory is truncated 
up to second order in spatial derivatives^. This action differs from the Arnowitt-Deser- 
Misner (ADM) formulation of Einstein's general relativity unless A = 1. Indeed, the 
presence of the constant A breaks the gauge symmetry of general space-time diffeomor- 
phisms and due to this it was believed that the effective theory (jl.ip reproduces general 
relativity only when A = 1. However, our analysis (7J showed that this is not the case; 
the dynamics of the action (11. ip is identical to the one of general relativity, regardless 
of the value of A. 

Let us take again and briefly the main point of the approach of Ref. [7J. When 
the preservation in time of the Hamiltonian constraint, which becomes a second-class 
constraint due to the reduced gauge symmetry group of Hofava models, is imposed, a 
new constraint emerges, 



Throughout this paper we shall call constraint C to the condition emerging when the 
preservation of the Hamiltonian constraint is imposed. Here 7r is the trace of the mo- 
mentum conjugated to gij. From Eq. (II. 2p and the asymptotic behavior of the field 

1 There have been also several analyses on exact classical solutions, some of them are [8]. 

2 The Hamiltonian analysis for the projectable counterpart of the action (jl.lj) was done in Ref. [5]. 
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variables it is straightforward to derive a global condition, which is 



/ 



d 3 x 



N 2 (Vi7cf = 0. 



(1.3) 



This equation leads unavoidably to the condition it = at each point of space-time. In 
this point we would like to remark that the possibility of iV = is discarded at the very 
definition of the action (II. ip . whose kinetic term depends on N as ~ 1/N. The same 
restriction must be imposed on the Hamiltonian, otherwise it lacks its relation to the 
action f ll.ip . Another important consideration is that this result is independent of the 
role one assigns to the variable iV in the canonical analysis: the result is the same if N is 
taken as a canonical variable or as a Lagrange multiplier, as we are going to see in Sec. 
II. Therefore, constraint C becomes completely equivalent to 7r = in the A-R-model, 
such that the kinetic trace-term drops out from the action (11.11) . We stress that n = 
is not an additional constraint put by hand. Its preservation in time yields an elliptic 
partial differential equation (PDE) for N, such that Dirac's algorithm stops here. 

When one focuses the problem from the side of general relativity, it turns out that 
7i = is a partial gauge fixing condition that can be imposed to fix the part of the 
gauge symmetry that mixes time with space and that is absent in the Hofava theory. 
When this partial gauge fixing condition is imposed on the side of general relativity, it 
becomes clear the equivalence between it and the A.R-model. 

Naturally, one question arises: can the same approach of the A-R-model be extended 
to other models with higher-derivative potentials?. The aim of this paper is to give an 
answer to this question by performing completely Dirac's algorithm of constraints for 
a model with a potential that goes one step beyond the previous A-R-model: we shall 
consider a potential with a R 2 term, which we shall call the A.R 2 -model. Thus, we shall 
study a fourth-order potential having a mathematical structure that is closer to the full 
Hofava theory than the A-R-model. 

We give an advance of our main results: Dirac's program can be performed success- 
fully for the A.R 2 -model, showing in this way the closure of the algebra. Constraint C 
is again solved completely as an elliptic PDE for tt, but now the solutions depend in 
general on N. This feature departs from the A-R-model and leads to the presence of one 
extra scalar degree, which is odd in the number of initial data it propagates. We also 
show a perturbative scheme under which this extra mode decouples at low energies. 

In Ref. jl] an extra, odd scalar mode was also found for a model with a RijR^ term 
in the potential. By making a perturbative analysis, those authors argued that the extra 
mode makes the original Hofava theory inconsistent. In particular, they concluded that 
the extra mode presents very fast instabilities and strong coupling at the IR. However, 
from their analysis one cannot arrive at that conclusion since their specific studies about 
the propagation of the extra mode were restricted to the A-R-model. Indeed, the main 
analysis of Ref. jl] concerning the perturbative propagation of the extra mode was based 
on the Stiickelberg formalism, which allows to cast the original Hofava model as a totally 
covariant theory using an extra scalar field. For the perturbative analysis those authors 
only took the reduced action given in Eq. (32) of their paper, which does not consider 
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the dynamics of the graviton because the 4D scalar curvature term is not taken into 
account. When one adds the 4D scalar curvature term and imposes the unitary gauge it 
turns out that the resulting action is exactly the A-R-model. Therefore, the perturbative 
analysis done in Ref. [I] cannot yield any information about the extra mode since it 
was performed at the order of the A-R-model and this theory does not possesses f£j|. 
Of course, at the moment those authors presented their study it seemed a consistent 
and convenient approach to analyze the extra mode using only the effective A-R-model 
since the presence of the constant A apparently entails a deviation from general relativity. 
After [7] we know that there is not an extra mode nor a deviation from general relativity 
for any value of A. 

This paper is organized as follows: in Sec. II we perform completely Dirac's algorithm 
for the A.R 2 -model. In Sec. Ill we carry out a perturbative analysis for the constraint C 
and its preservation in time. To achieve this we take the approach of the weakly varying 
background used in Ref. [I], showing the decoupling of the extra mode at linear oder 
in perturbation and zeroth- and first-order in derivatives of the background. In Sec. IV 
we present our conclusions. 



2 Analysis of the constraints 

The action of the A.R 2 -model, in ADM variables, is 

S = J dbd z x^N{G m K i:j K kl + R + (R 2 ) , (2.1) 

where -R is the spatial Ricci scalar, ( is a coupling constant and the kinetic variable 
is the extrinsic curvature defined by 

K* = ^(fti - 2V (4 JVfl) . (2.2) 
Of course, this requires that N ^ in all space-time. G^ is a four-indices metric, 

G w = }_ {g * kgJ i + gi i gjk) _ Xgijgk i (2 3) 

For the case of A ^ 1/3 the inverse of G % i kl is given by 

Gijki = ^(9ik9ji + 9u9jk) - 3^ _ ^ 9ij9ki ■ (2.4) 



3 In addition, in Ref. [4] a preliminary perturbative analysis was done besides the Stiickelberg ap- 
proach. There those authors presented a dispersion relation for the extra mode. However, that study 
was again restricted to the Ai?-model, as they indicated in their paper. Indeed, If one takes Eqs. (21) 
and (22) of Ref. [1] and for the background imposes K = 0, which is the only solution of the C constraint 
of the Ai?-model, one finds from Eq. (22) that k = and from Eq. (21) n = 0, k being the perturbation 
of K and n the one of N. This result is consistent with the fact that the AJ?-model has no an extra 
mode. Hence, the dispersion relation (23) and all the analysis that follows is invalid. 
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The momenta conjugated to N and Ni vanish automatically for the Lagrangian given 
in Eq. ( 12. ip . The same happens in general relativity and in any Hofava model. Because 
of this N and Ni are usually considered as Lagrange multipliers, as we did for the XR- 
model in Ref. [TJ. However, for the A.R 2 -model we find it convenient to keep the lapse 
function N, together with its canonical conjugated, as part of the canonical variables, 
whereas the shift functions Ni are regarded as Lagrange multipliers and their conjugated 
momenta are dropped out from the canonical space. We denote by <fi the momentum 
conjugated to N and by tt^ the one conjugated to gfy, which is found to be 

^ = G^ kl K kl . (2.5) 

This formula is the same for any theory having the same kinetic term of the Hofava 
theory. The velocities can be solved explicitly in terms of the momentum 7r u whenever 
A / 1/3. 

We assume that the whole spatial manifold is noncompact. All the configurations are 
fixed at infinity and correspond to Minkowski space-time. We also assume that the field 
variables have the same asymptotic behavior of general relativity: in asymptotically flat 
coordinates, as r approaches to infinity, they behave as 

g .. = 5 .. + 0(r -i) ; 7r« = 0(r- 2 ), 

(2.6) 

N = l + 0(r~ l ), Ni = 0{r~ l ). 

We need to clarify that here the notation 0(r~ a ) means that the variable asymptotically 
receives contributions at most of order r~ a , but they can be effectively lower. For 
example, in the XR- model if 7r|oo = then it is zero everywhere, hence it does not 
receive any contribution. 

By performing the Legendre transformation, assuming the A ^ 1/3 condition, we 
obtain the Hamiltonian 

H = (NH + NiW + acj)) , 

W EE -2V^\ 

where we use the brackets ( ) to denote integration over all the spatial manifold, (J 7 ) = 
J J-d 3 x, and the standard definition for covariant derivatives on densities is understood. 
We have added up the primary constraint 0, a being its Lagrange multiplier, such that 
the above is the total Hamiltonian. Under spatial coordinate transformations <ft behaves 
as a scalar-density, hence a is a scalar. Notice that this Hamiltonian has been deduced 
under the assumption N ^ 0, thus we must keep this restriction in the Hamiltonian in 
order to be able to associate it to the A_R 2 -model. 

Now we start the Dirac's algorithm to extract all the constraints of the theory and to 
ensure their preservation in time. The condition TV = emerges as a primary constraint 
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when one takes variations of the canonical action w.r.t. iVj. The preservation in time 
of this constraint is ensured due to its role as the generator of the spatial coordinate 
transformations [10] . Hence, as the next step we must demand the preservation in time 
of the primary constraint <f), which yields 

{^H} = -H. (2.8) 

Therefore H is a secondary constraint of the theory. 

The preservation in time of % requires to compute only the bracket of H with 
itself, {H, (N%)}. We obtain that the vanishing of this bracket leads to the secondary 
constraint C = 0, where 

C = Vi(iVV), (2.9) 
7 i = /VV^vr/p) + 2Cn ij V 3 R (2.10) 

and in addition n = g^n^ , p = 1 + 2(,R and = (A — 1)/(3A — 1). Notice that 
expression (12.101) is regular in the limit p = 0. Constraint C is the extension of the extra 
constraint found in the \R- model [U [7J. It arises as a consequence of the lacking of the 
general space-time diffeomorphisms as gauge symmetries, which leads to the second-class 
behavior of the Hamiltonian constraint %. 

Having found the constraint C, in order to obtain a consistent closure of the algebra 
it becomes crucial to determine the appropriate variable for which its general solution 
can be found in a closed way. We recall that we have adopted a scheme under which 
N, as well as gij and ir %3 , are canonical variables, hence all of them are treated on the 
same footing. This allows us to solve, in principle, any of these variables in terms of the 
other ones. 

The first thing to do is contrasting with the Ai?-model: if the vector-density field 
7* has the form 7* = % 2 VV f° r some nonzero function x an d a function ip vanishing 
at infinity, then, by following the same steps done in Ref. [7J, one easily obtains that 
the constraint C does not give any information about N, it instead implies p = 0. This 
is the case for the A-R-model, where 7* = /3VV. For the Ai? 2 -model, the first term in 
7* (I2.10p has exactly the required form, but the last term does not, hence we are faced 
with an obstruction to apply this procedure. This leads us to conclude that there are 
solutions of C with nonvanishing vector-density field 7* which are iV-dependent. 

The next question we might ask is if the general solution of C can be cast as a 
condition for N. We are going to show that the answer is negative by making use of 
the asymptotic expansion of the constraint C for large r. This is based on the fact 
that the values of the field variables at infinity are fixed and correspond exactly to 
Minkowski space-time. Thus, the degrees of freedom of any field variable are expressed 
asymptotically in the coefficient of its first corrections. By using the expansion (12.61) 
in (12.101) we obtain that 7* = 0(r~ 3 ) asymptotically. Then, it is clear that the leading 
term in the asymptotic expansion of C, which is of order 0(r~ A ), comes from the pure 
divergence of 7* and does not involve the variable N. Consequently, the general solution 
of constraint C necessarily restricts 7*, hence and 7r y , at least asymptotically. 
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In order to solve constraint C in the most general way, we cast it as an equation for 
7r. It results that for this variable C is an elliptic PDE compatible with the asymptotic 
conditions. To see this it is convenient to make the covariant decomposition of 7r u 
between its trace and traceless parts, 

71^ = 7^' + ^%. (2.11) 

Using this decomposition, the constraint C given in Eqs. (I2.9p and ( 12 . 1 j) takes the form 

[/3iVy V 2 + (0{N 2 p 2 ) + \N 2 d l p 2 ) V, + |V,(iV 2 9V 2 )] (tt/p) = -2CV,(iV 2 7r^) , 

(2-12) 

which, for A ^ 1, is an elliptic PDE for n/ p with an iV-dependent source term. There 
always exists a unique solution for it satisfying the boundary conditions. The interpre- 
tation of C we make here is similar to its role in the Ai?-model, with the difference that 
here the solution for 7r depends on N. This is the main reason why it is convenient to 
regard N as a canonical variable in the Ai? 2 -model. 

Now we impose the preservation in time of the constraint C (12.91) . which is established 
by the condition {C, (N7i + &<j>)} = 0. We evaluate this condition by direct computa- 
tions and simplifying the result using all the known constraints. This yields an equation 
for a: 

V'Viia/N) = -fVi [iV 2 (/3p 2 ^ + (& + ACg-^p^)] , (2.13) 

where 

V = ivv/vs, 

B l = [N(R + (R 2 ) - V 2 (Np)} V l R + 2 [V*V J '(Np) - NpW] VjR 
+V l [NR(R + (R 2 )} + f3p 2 V i (NR 2 /p) 

+ [ 2 (7r w V fc V,iV - (5V 2 {Nrx) - Ng klmn n M R mn ) V^vr/p) 
9 

+ P V i (N7rg iklm ^ k R lm /p 2 ) - pW^ix^WkWiN/p 2 ) + f3pV l (rrV 2 (Nrr)/p 2 )] , 
V ij = n ij (Nir kl R H - 7t kl V k ViN + /3V 2 (Ntt)) - NRn k ^TT j \ . 

(2.14) 

Equation (I2.13P defines a first-order PDE for a in any open set in which the vector 
field V 1 is nonvanishing. The standard theorems about PDEs ensure the existence and 
uniqueness of the solution once a suitable boundary condition has been provided. Since 
l^Vj is a flow-operator, the boundary data compatible with it consists of giving the 
value of a/N on a noncharacteristic 2D surface, that is, a surface that is nonparallel to 
V % at each point. Notice that, due to constraint C, the vector field V 1 is divergenceless, 
and by the required boundary conditions it vanishes at infinity with a decay of order 
0(r~ 3 ). As a consequence, it can be shown that its flow cannot reach infinity. Further 
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details about the flow of V % and the way in which the Eq. ( 12.130 must be solved can be 
found in the Appendix. 

We wish to clarify the compatibility between the generic solutions of Eq. ( 12.130 and 
the expected asymptotic behavior of the field variables. To this end we first point out 
that the canonical equation of motion 8S/8<fi yields 

N = a, (2.15) 

which gives the evolution of the variable ./V once a has been solved from the Eq. (12.131) . 
From Eqs. (12.60 and (12.150 we deduce that asymptotically a behave as a = 0{r~ l ), 
such that the l.h.s. of Eq. (12.130 . which involves a, has an asymptotic decay of order 
0(r~ 5 ). In the r.h.s., the leading term comes from the first term of A\ which also leads 
to a decay of order C(r~ 5 ) in the r.h.s. of Eq. (12.130 . Thus we see that the Eq. (12.130 
imposes consistent restrictions on the multiplier a at any order asymptotically. We 
conclude that Eq. (12.130 can be consistently used to solve for the Lagrange multiplier 
a, ending Dirac's algorithm with this step. 

We have then ended up with the first-class constraint W = (3 constraints) and 
with the second-class ones <p — 0, 7-L = and C = (3 constraints). The canonical 
variables are {N, <fr, gij, 7r u } (14 variables). The number of physical degrees of freedom 
in the canonical space is given by 

(# Physical D.O.F.) = 

V ; (2.16) 

(# Canonical var.) - 2x (# 1 st class const.) — (# 2 nd class const.) = 5 . 

Four of these 5 degrees of freedom correspond to the propagation of the graviton in 
the phase space. The remaining degree is an odd scalar mode, that is, a scalar field 
that propagates only one Cauchy field as initial data, but not its time-derivative. This 
extra mode is represented by the canonical variable ./V and the corresponding evolution 
equation is given in (12.150 . The general analysis about the propagation of the extra 
mode is as follows: once a has been solved from Eq. (12.130 on a t = constant slice 
A^3, then the evolution of iV is determined by Eq. (12.150 together with the initial data 
N(t = 0, x) provided. This initial data represents extra information in the sense that it 
does not arise in general relativity. 

It is worth mentioning that our analysis is based on the fact that the solution for the 
multiplier a is constructed by following the integral lines of the vector field V i . There 
could be, of course, some configurations for which V % vanishes in all points of space- 
time. Those points of the phase space are very special since they solve automatically 
the constraint C for arbitrary N and also the Eq. (12.130 cannot be used to solve for 
a in any open set. As an example, this is the case of all the configurations satisfying 
R = — 1/2( (p = 0) at each point. The algebra of constraints on those configurations 
deserves a detailed analysis besides the study we present here. 

As a check of consistency we now go to the limit ( = and compare with the 
A-R-model. In this limit the Hamiltonian constraint of the Ai?-model is recovered from 
(12.70 and Y given in Eq. (12.100 is reduced to Y = /3VV, such that C also takes the 
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same form of the A-R-model and 7r = everywhere is its only solution satisfying (12. 6p 
asymptotically. However, note that in the scheme we have adopted in this paper N is 
a canonical variable, as a consequence Eq. (12.131) must still be analyzed. This equation 
simplifies greatly since Y vanishes in the ( = limit due to the C constraint. The 
l.h.s. of Eq. ( I2.13P vanishes, consequently we do not get from it an equation for ex, but 
an extra constraint, 

0V, [iV 2 V (V 2 - R) N] = . (2.17) 

This constraint is solved in the same way one solves C in the A-R-model: the scalar (V 2 — 
R)N is of order 0(r~ 3 ) asymptotically; one multiplies Eq. (I2.17P by it and integrates 
over all the spatial manifold. From the resulting integral, after an integration by parts, 
one deduces that the only solution is (for A ^ 1) 

(V 2 -R)N = (2.18) 

at each space-time point, which is exactly the equation that determines N in the A-R- 
model [7]. In that reference this equation was obtained by requiring the preservation 
in time of the constraint 7r = 0. In the ( = limit of the A-R 2 -model the mechanism 
works analogously, but being careful in identifying the Eq. (12.131) as an equation for N, 
while a decouples completely from it. Equation (I2.18P must be regarded as a further 
extra constraint of the theory, such that the number of degrees of freedom is reduced in 
one. In this scheme the Lagrange multiplier a will be determined once the preservation 
in time of the constraint (12.181) is imposed (or directly by Eq. (I2.15p . which yields the 
same result). 



3 Perturbative analysis on a weakly varying back- 
ground 

The general analysis done in the previous section has revealed the presence of an ex- 
tra, odd, scalar mode. However, this extra mode is not found when one analyses the 
dynamics of the lowest-order effective theory, the A-R-model, taken as an autonomous 
theory [7|. Therefore, it becomes crucial to determine if the extra mode we have found 
in the A_R 2 -model decouples at low energies, such that the physics of the A-R-model is 
smoothly recovered in the IR limit. 

A perturbative analysis is a good approach to elucidate this point. The conditions 
that play a role in the coupling/ decoupling of the extra mode are the constraint C and 
its preservation in time, which yields the Eq. (12.131) . We may concentrate ourselves only 
in the perturbative expansions of these conditions and in particular we must study the 
perturbative behavior of the variables N and a. With regard to the background, our 
main interest is not in the specific evolution of the modes around a given vacuum, we 
are more interested in seeing the IR decoupling of the extra mode for a broad class 
of backgrounds. To achieve this, we adopt the scheme used in Ref. jlj: instead of 
dealing with a specific background, we assume that a solution of the constraints and the 
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equations of motion, denoted by {g^-, , N, a}, exists and expand the fields around the 
solution up to linear order, 

9ij = 9ij + hij , ir ij = n ij + //' , 

N = N + n, a = a + r. 

It is supposed that the solution depends both on space and time, hence all of the fields 
of the background are in principle different from zero. In order to further simplify the 
analysis, it is also assumed that the background is a weakly varying configuration, in 
particular when it is compared with the space-time scales in which the perturbations 
varies. Specifically, the background metric, together with N, changes at a typical space- 
time scale L. This gives ~ L' 1 , R ~ L~ 2 , diN ~ L~ x and a ~ L^ 1 . Under this 
assumption one can assign a weight to each term in the perturbative expansions of the 
constraints/equations of motion according to its order in L^ 1 . 

We perform a linear-order perturbation on the constraint C given in Eqs. (12.91) and 
fl2TTUD and obtain 

C = d i (N 2 Y)+2d i (N?n) , (3.2) 

where 

ryi _ ryi(0) _|_ ryi(l) _|_ ryi(2) _|_ y(3) _|_ ^l(4) ^ 

y(o) = 

7 i(i) = pTt^hk! - 2/3CTrd i 5R + 2C7t ij dj5R - /37tg ij 8r jk k , 

7^ 2 ) = 2/3(RV i p + 20CV i it5R - pVjirh ij + PV^ kl h kl , 

(3.3) 

7 i(s) = 2 ( {f3R=n kl V i h kl - (5d l Rp + p ij djR - 0R7rg ij 5T jk k ) , 

7 *(4) = 2 ( (pRV^hki - 0m kl h H - PRVjTrtii + 0TtdjRh ij ) , 

5T i3 k = V (i ^) fc -iV%, 

= V, V ; //'' - V 2 /i - . 

In these expressions all spatial indices are raised/lowered with cjij, h = g^hij and p = 
(jijP 13 ' . In (13.31) we have grouped^ the terms of 7* according to their order in L _1 and all 
the terms that are of order zero in the perturbative variables have been omitted because 
the background is a solution of the constraints of the theory. 

In the expansion (13. 2p the perturbative variable n arises only in the second term. 
Thus, the term of lowest order in L _1 containing n is +2/3iVV*7f<9jn, which is of order 

4 Actually, the classification is a bit disordered since covariant derivatives of the perturbative variables 
are not of homogeneous order in L _1 : terms having the background connection are of one order higher 
than the pure derivative term. This only means that in the expansion of 7* these terms do not belong 
to the order they have been written, but to the next order. 
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L~ 2 . Therefore, in the linear-order perturbative analysis the variable n decouples from 
the constraint C if we expand it up to order L^ 1 , resulting in a constraint for the variables 

hij and p % \ 

C = C<°) + C (1) , (3.4) 

where 

C<® = ^N^id'p, (3.5) 
C« = ^jV 2 ^ln(iV7v^)^P 

+N 2 (pjPdi&hia - 2P(fid i d t 5R + 2(7t ij d i d j 5R - >r,ir<SV, f ') . (3.6) 

We notice that, when expanded up to order zero in L -1 , constraint C acquires exactly 
the same perturbative form of the A-R-model: it becomes = and its only solution 
vanishing at infinity is p = everywhere. 

Now we move to the Eq. ( 12.1 3ft . For our purposes it is enough to expand this equation 
up to linear order in L _1 . This yield^] 

&diV*n + di \n{^JgN 2 )d i 'd j d j n + 2(Nd i d i V 2 5R + 2((6diN + Nd, t In ^d'&d^R 

-N&dtfR - (49W + Ndi In y/g)di5R - d^g^did^r^ 

3Ntt a 



2 "' V [^{d'difp + ^did^dup] 



2(3A - 1)5 yfiN 



d l dip = 0. 
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(3.7) 

As we may see, the perturbative multiplier r does not arise in this equation, which is 
a fourth-order elliptic PDE for n. This result is consistent with the decoupling of the 
extra mode from the constraint C. Equation (13.71) reduces in one the number of physical 
degrees in the phase space, showing the consistent decoupling of the extra mode at linear 
order in perturbation theory and zeroth and linear order in L~ l . 

Let us discuss the smooth decoupling of the extra mode from the point of view of the 
initial Cauchy data. As we pointed out in previous section, the initial data N(t = 0, x) 
is needed in the general analysis for the complete evolution of the system. However, in 
perturbation theory the Lagrange multiplier r is fixed by Eq. ( 12. 13f) only at second order 
in L _1 and higher, whereas at zeroth and first order the Eq. (12.131) is an elliptic PDE 
for iV that fixes it for any time. This implies that at zeroth and first order in L _1 the 
equation of motion (12 . 1 5[) must not be regarded as an independent evolution equation 
since the evolution of N is derived from the other canonical variables. Equation (12.151) is 
instead the equation that determines r at these orders. Therefore, the evolution of the 
system at zeroth and first order in L~ x does not require N(t = 0, x) as an independent 
initial data. It is involved in the initial value problem from second order in L~ l and 
higher. 



5 The expansion has been computed assuming A ^ 1. There is not any contribution to the Eq. (| 2 . 1 3 1) 
at the orders we are considering if A = 1. 
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4 Conclusions 



We have carried out the Hamiltonian analysis for a Hofava model with a potential 
composed of R- and i? 2 -terms. We have performed Dirac's algorithm for extracting all 
the constraint of the theory and ensuring their preservation in time, showing in this way 
the closure of the algebra of constraints. 

The constraint emerging when the preservation in time of the Hamiltonian is imposed 
can be solved for the variable 7r = g^ir 1 -* in a closed way, since this constraint is an elliptic 
PDE for 7r. Then, the preservation in time of this constraint is imposed and it turns out 
that Dirac's algorithm ends at this step since this yields a first-order PDE for a Lagrange 
multiplier of the theory. This equation is a flow equation along the integral lines of a 
divergenceless vector field and its solutions are compatible with the expected asymptotic 
behavior of all the fields. Our analysis provides non trivial evidence about the consistent 
structure of constraints for Hofava models that go beyond the lowest-order potential. 

We have seen that the inclusion of the R 2 term in the potential gives rise to an 
additional scalar degree of freedom that is absent in the XR- model analyzed in Ref. [TJ. 
This scalar is an odd mode in the sense that it propagates only one field as Cauchy data. 
This agrees with the extra mode found in Refs. [H H] . We have also performed a linear- 
order perturbative analysis on a weakly varying background, showing the decoupling of 
the extra mode for zeroth- and first-derivatives of the background. We have discussed 
how the initial data associated to the extra mode smoothly decouples from the initial 
value problem of the whole system. 

The physical propagation of the odd extra mode is an analysis that must be done 
carefully, due to its peculiarity of having an evolution equation of first order in time- 
derivative, as was also noticed in Ref. [I]. In any case, the analysis should effectively 
include higher-order terms in the potential. 

Our results indicate the classical consistency of the truncated Ai? 2 -model. However, 
the validity of the full nonprojectable Hofava theory still remains unanswered. In fact, 
as we have mentioned, even in perturbation theory the propagation of the extra mode 
has not been completely analyzed. Moreover, the quantization of the Hofava theory is 
a challenging task due to the presence of very involved second-class constraints. 

Appendix: Analysis of the equation for a 

The main difference between the Eq. ( 12.1 3p for a and the standard constraints of general 
relativity is that Eq. (12. 13[) is not an elliptic PDE. This deserves some discussion about 
the nature of the differential operator of Eq. ( 12 . 13[) and the boundary data compatible 
with it. 

We consider the existence and uniqueness of the solutions of Eq. (I2.13P on a t = 
constant slice Ai^, which is assumed to be a Riemannian manifold topologically equiv- 
alent to M 3 . We may rewrite (I2.13P as 

V%(a/N) = J , (A.l) 
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where the source J is independent of a. All fields are assumed to be smooth geometrical 
objects on 

The l.h.s. of Eq. flA.lj) is the application of the divergenceless (constraint C) vector 
field V, whose components are V 1 = N 2 Y/y/g, on the scalar field a/N. Therefore, 
in order to analyze the Eq. ( lA.ip it is important to describe the flow generated by V. 



To this end we shall consider some mild assumptions on V, which in particular may 
be imposed on the initial data over the slice t = when formulating the initial value 
problem to the field equations of the full theory: We assume that the zeros of V belong 
to a 2D surface £ C A^3, where £ splits A4 3 into two open sets denoted by M. + and 
A4-. £ may be extended to infinity or it may be a compact smooth surface without 
boundaries. 

The flow determined by the vector field V is obtained as the solution of the ODEs 

— = V(*-00), (A.2) 
as 

where x % are local coordinates on A4 3 . The solution x % = x*(s,Si,s 2 ) depends on the 
flow parameter s and si, s%, where x*(0, S2, S2) is the local parametric description of S. 
In turn, the Eq. flA.lj) can be cast as 



d(a/N) , , k , 

K ' J =J(x'(s, Sl ,S2 ■ A.3 
as 



The system of ODEs (I A. 21) - (1A.3j) is equivalent to the PDE (lA.lj) in a neighborhood of 



£. The change of variables x % -H- (s, si, S2) is a diffeomorphism in this neighborhood. 
Let us denote by u the unit vector field orthogonal to S and oriented inwardly A4+. 

Let 

[/± = {pGE : and U = {p € E : F ■ M = 0} . (A.4) 

The continuity assumption ensures that any curve on E from pi G U + to P2 G nec- 
essarily intersects Uq. Notice that U + by definition is automatically a noncharacteristic 
surface with respect to the vector field V . Our aim is to show that U + is the appro- 
priated surface to give the boundary condition on a/N in order to solve uniquely the 
Eq. ( TA~Tj) on M+. 

There are two important facts resulting as consequences of the constraint on V and 
the settings we are considering. The first one concerns to what extent the flow of V 
can reach the infinity. We remark that V is divergenceless and its asymptotic behavior 
(0(r~ 3 )) ensures that its flux at infinity is always zero. Thus, one can take a Gaussian 
surface in A4+ with the form of a tube parallel to the flow lines of V and having one of 
its caps coinciding with an arbitrary subset Q C U+ of nonzero measure and the other 
cap at infinity. By considering the flux of V on the whole Gaussian surface one obtains 
automatically that the flux on Q vanishes. Therefore, it turns out that if we take points 
on U + , only a set of measure zero, if any, may have associated curves x l (s, si, §2) ending 
at infinity. The main flow of V always ends at LL or Uq [see Fig. [TJa]. 
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Figure 1: (a) A typical flow of V. (b) If a zero of V was excluded from E then in the 
neighborhood of the zero there could be a set of loops disconnected to E. 



The second point is the fact that the flow of V starting on points of U+ defines a 
map x l (s, s\, 82) '■ U + — > M. + which is a surjective map on M. + . In this point we remark 
that, since E contains all the zeros of V, if a loop 011 M3 circulates around a zero we 
avoid situations like the one showed in Fig. [T] (b). In that case the flow starting at U+ 
does not span completely A4 + , hence the flow does not determine a/N inside the loop. 
In our case any loop on jW 3 circulating a zero necessarily intersects E. 

After these considerations we have arrived to our main conclusion about the equa- 
tion for the Lagrange multiplier a. The standard theorems on ODEs ensure the local 
existence and uniqueness of the system ( 1A.2I) - ( 1A.3j) and hence of the PDE ( 1A.1I) on 



given an initial data / + on any open set Q G U+, there always exists a unique 
smooth solution a/N of Eq. (lA.lj) in a neighborhood Q satisfying (a/N)\a = /+. 



Our qualitative analysis suggests that the solution may be extended globally. How- 
ever, in order to address that problem a rigorous discussion of the propagation around 
the zeros of the vector field V must be considered. We will analyze this problem else- 
where. 

The Dirac algorithm for the Ai? 2 -model ends at this step since the Eq. ( 12 . 1 3[) is 
solved for the Lagrange multiplier a. However, it is interesting to notice that a is 
not completely determined by Eq. fl 2 . 1 3 p : the initial data / + at U + is arbitrary and 
the solution a depends continuously on it. This means that the associated primary 
constraint is a mixture of first- and second-class constraints. The gauge symmetry 
generated by the first-class part has a gauge parameter depending only on the local 
coordinates on E. In particular, this lower-dimensional gauge symmetry may be used to 
fix a boundary condition on E for the scalar field N. Besides this subtlety, N describes 
a half-degree of freedom that is absent in general relativity. 
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